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The Kondo system CeRu 4 Sn 6 shows a strong anisotropy in its electric, optic and magnetic prop¬ 
erties. We employ density functional theory plus dynamical mean field theory and show that the 
predominant Ce-f state has total angular moment J — and z-component mj — ±1/2 in agree¬ 
ment with recent X-ray absorption experiments. Even though CeRu 4 Sn 6 has the direct gap of a 
Kondo insulator through most of the Brillouin zone it remains weakly metallic. This is because of 
(i) a band crossing in the ^-direction and (ii) a negative indirect gap. 


PACS numbers: 71.27.+a, 31.15.V- 

In /-electron systems such as CeRu 4 Sn 6 we have 
the generic situation that narrow, strongly interacting 
/ bands hybridize with weakly correlated conduction 
bands, as exemplified by the periodic Anderson model. 
This gives rise to a hybridization gap at the crossing 
of / and conduction band, already without /-/ interac¬ 
tion e.g. in density functional theory (DFT). Electronic 
correlations strongly renormalize this gap, resulting in 
a gapped Kondo resonance which can be understood in 
a quasiparticle picture. If the gap is at the Fermi level 
one has a Kondo insulator; otherwise one has a metal 
with a typically heavy effective mas^^^^. Even more com¬ 
plicated is the situation in anisotropic Kondo insulators 
such as CeNiSiP. Here, the physical properties such 
as the susceptibility and the conductivity are strongly 
anisotropic. One possible explanation is a hybridization 
gap with nodeP® so that in some directions the bands 
cross the Fermi level whereas there is a gap in other parts 
of the Brillouin zone. Besides, also a V-shaped DOP 
or extrinsic effects such as off stoichiometry, impuritie^ 
or topological surfaces states have been considered as a 
microscopic origin of the residual metallicity in some di¬ 
rections. 

Such an anisotropy has also been observed in 
CeRu 4 Sn 6 single crystalP with a sirmler, tetrago¬ 
nal I42m crystal struct urd^^^ (Fig. the opti¬ 

cal conductivity shows a weak Drude-like feature in 
the a — a plane whereas it has a dip at low fre¬ 
quencies in the c-direct iorP^. This has been con¬ 
firmed by density functional theo ry plus dynamical mean 
field theory (DFT±DMFT]P^^ calculations of t he op- 
tical conductivit}ff^. Likewise the thermopowei^^^^, 
resistivit}!^ and magnetic propertie^^ are strongly 
anisotropic. 

In the present paper, we discuss the DFT±DMFT cal¬ 
culations in more detail. In particular we show that 
electronic correlations and the Kondo effect reverse the 
crystal field splitting between J = 5/2 \mj\ = 1/2 and 
J = bl2 \mj\ =3/2 Ce-f states. Here, J is the total 
angular moment and mj its ^-component. Consequently 
the \mj\ = 1/2 states get predominantly occupied with 
decreasing temperature (increasing correlation effects). 



FIG. 1: Crystal structure of CeRu 4 Sn 6 . 


Nonetheless, the evolving Kondo resonance at the Fermi 
level also contains \mj\ = 3/2 states. This DFT±DMFT 
finding agrees with recent x-ray experiments^ which best 
fit to \mj\ = 1/2. We also discuss why CeRu 4 Sn 6 re¬ 
mains metallic. In section [I] we present the DFT band- 
structure and Wannier orbitals. In Section El the or- 
bitally resolved DMFT spectra and self energies and their 
temperature-dependence is discussed. Finally, Section [TTI| 
summarizes our results. 


I. DENSITY FUNCTIONAL THEORY 

Starting point is the experimental tetragonal lA2m 
structure of CeRu 4 Sn 6 with lattice parameters a = 
6.8810 A, and c = 9.7520 Here, the Ce atoms are 
surrounded by a ’’cage” with four nearest neighbor Ru 
atoms and four next nearest neighbor Sn atoms which 
are only 2% further away, see Fig. [I] Our first step is 
to calculate the DFT electronic structure with Wien2K!^ 
using the PBE-GGA potential and a fine /c-mesh of 10000 
points to obtain accurate results. Spin-orbit effects are 
important due to the heavy atoms and hence included. 














2 


The DFT partial densities of states for CeRu 4 Sn 6 are 
shown in Fig. The calculation reveals a mixed man¬ 
ifold around the Fermi level with main contributions of 
Ce-/, Ru-d and Sn-p characteJ^. The narrow /-bands 
are essentially placed between 0.1 and 0.7 eV above the 
Fermi level in DFT. Because of the small hybridization 
tail below the Fermi level, the DFT /-filling is nonethe¬ 
less surprisingly large, i.e., 0.7 electrons per Ce atom. 

Fig. 1^ (bottom) shows that the spin-orbit coupling sep¬ 
arates the Ce-/ states into J = 5/2 states, whose center 
of gravity is approximately around 0.2 eV, and J = 7/2 
states around 0.5 eV. The crystal field further splits the 
J = 5/2 manifold into mj = ±5/2 character at higher 
energy and mj = ±3/2 as well as mj = ±1/2 states at 
essentially the same energy. The crystal field splitting 
between the latter is only ^ 0.05 eV with mj = ±3/2 
below mj = ±1/2. This refiects in the small splitting 
in the /-level at ~ 0.2 eV in Fig. (bottom). Nonethe¬ 
less, the largest contribution below the Fermi level has 
(J = 5/2, mj = 1/2) character because these states hy¬ 
bridize more strongly with Ru-d than the other / orbitals 
[the integrated weight below the Fermi level is larger for 
the mj = ±1/2 orbital in Fig. (bottom)]. Apart from 
the / partial density of states, the major manifold in 
Fig-i is Ru-d with traces of Sn-p. 

The corresponding bandstructure of CeRu 4 Sn 6 is 
shown in Fig There is a small dip above the Fermi 
level. This is due to a DFT direct gap between the 
highest lying occupied band (mainly Ru-d) and the low¬ 
est (unoccupied) / states with mixed mj = ±3/2 and 
mj = ±1/2 character. As one can see in Fig[^ the va¬ 
lence and conduction band both touch the Fermi level 
however so that the indirect DFT band gap is zero. Also 
note that the lack of inversion symmetry around the Ce 
atoms lifts the degeneracy between the ±m j components. 

Next, we do a projection onto maximally localized 
Wannier orbitalP^ using the Wien2Wanniep^ package. 
In this material, the spin-orbit coupling is crucial which 
makes the identification of adequate Wannier orbitals for 
modelling the low energy degrees of freedom more dif¬ 
ficult. The decisive aspect is the choice of the atomic 
orbitals and for the Wien2Wannier localization to 
start with. For example, for (J = bl2^mj = ±3/2) we 
start with 



where l^s+i) 1 ^ 3 + 2 ) denote the 

\L = 3,771^ = ±1,5' = 1/2,rUg = ±1/2) and 

|L = 3, = ±2, 5 = —1/2, rus = ±1/2) in the (Im)- 

basis, respectively; the prefactors are the Clebsch- 
Gordan coefficients. Besides the Ce-/ orbitals, we 
included all 40 Ru-4d and 36 Sn-5p orbitals in the 
Wannierization. Fig. shows the most important 
Ce-J = 5/2 orbitals and their (fat-band) contribution 
to the bandstructur. The spread of the Ce-/ Wan¬ 
nier orbitals is: Os ,3 = 0.97 A, Osii = 1.10 A, 

9 5^ 9 9 5^ 9 



FIG. 2: Top: DFT partial density of states resolving the 
most important atom-wise contributions: Ce-/, Ru-d and Sn- 
p. Bottom: Zoom in and further differentiation of the f-states 
in the (J, mj)-basis. 

Cs ,5 = 1.83 At Ci ,7 = 1.93 At Qi .5 = 0.86 At 

Or I 3 = 1.32 A^, and Or ,1 = 0.99 A^. That of the 

Ru-d Wannier orbitals ranges from 2.32 to 5.79 A^, and 
that of the Sn-p is even larger with O from 5.89 to 18.14 

k\ 


II. DFT±DMFT 

We supplement the 90 x 90 Wannier Hamiltonian by 
a local interaction U = U' = 5.5 eV between all Ce- 
J = 5/2 states and subtract the double counting in the 
fully localized limifP^. The corresponding Hamiltonian is 
then solved by DMFT with Hirsch-Fye quantum Monte 
CarlcP^ simulations. Note that the Ce-J = 7/2 orbitals 
remain unoccupied well above the Fermi level so that we 
do not need to include them as interacting in DMFT. For 
the sake of simplicity we neglect Hand’s exchange which 
is justified since Ce only has about one /-electron and 
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employ the same DMFT self energy for ±mj. 

Let us start our discussion by looking at the orbitally- 
resolved local electronic densities. As already mentioned, 
in DFT the \mj\ = 1/2 orbital is slighly more occupied. 
Specifically the electron density of the Ce-J = 5/2 man¬ 
ifold is 

noFT = (^|mjl = l/2:^|mjl=3/25^|mjl=5/2) 

= (0.34 0.20 0.28) (2) 

( 3 ) 

Electronic correlations lead with decreasing temperature 
to higher \mj\ = 1/2 occupations whereas the \mj\ =5/2 
states get depopulated in DMFT: 


FIG. 3: DFT bandstructure around the Fermi level 
(solid,gray) perfectly reproduced by the Wannier orbital 
bands (dashed,green). Here and in the following figures, 
G = [0 0 0], A = [0 0 |], A1 = [1 0 |], A = [0 I 0], 
A2=[l lO], A3=[| I |]. 



(J = 5/2,mj = -h3/2) 



(J = 5/2,mj = -h5/2) 



FIG. 4: Wannier orbitals for the Ge-f manifold in GeRu 4 Sn 6 
and corresponding fat-band plots. For the charge density 
plots of the Wannier orbitals, the surface shows \w{r)\‘^ = 5 
Bohr“^ and the color show the sign of the phase 0, i.e. 
red (blue) for cos(0(re(r))) positive (negative). 


UDMFTiT =1160K) = (0.44 0.38 0.15), (4) 
UDMFTiT = 290 K) = (0.60 0.27 0.08). (5) 

The higher \mj\ = 1/2 occupation can be traced back 
to a reversal of the order of the crystal field levels. As 
discussed above, in DFT \mj\ = 3/2 is the lowest crystal 
field state, 0.05 eV below \mj\ = 1/2. To this differ¬ 
ence, we need now to add however the real part of the 
self energy, which altogether yields the so-called effective 
crystal field. The DMFT values (subtracting ja) for the 
three mj states are, also see Fig. 

5?E(T = 1160A:)-/i = (-0.61 - 0.26 0.91) eV, (6) 
5?E(T = 290A:)-/i = (-0.13 0.13 1.09) eV. (7) 

That is the \mj\ = 1/2 is shifted down and eventu¬ 
ally becomes the lowest effective crystal field state in 
DFT+DMFT. We attribute this to the larger hybridiza¬ 
tion between \mj\ = 1/2 and conduction electrons as 
well as to the Kondo effect (spin-flipping mj = ±1/2 
to mj = ±1/2 is possible for conduction electrons with 
angular momentum 1/2). 

The Kondo effect can be identified by the development 
of a Kondo resonance at the Fermi level when decreasing 
the temperature from T = 1160 K to 290 K, see Fig. 
Fig.Q shows a zoom in and an orbital breakup. Clearly, 
the Kondo resonance at the Fermi level is made up from 
mj = ±1/2 and mj = ±3/2 states; the mj = ±5/2 
states are at a somewhat higher energy due to the (cor¬ 
relation enhanced) crystal field splitting. 

For CeRu 4 Sn 6 , there has been a long debate over the 
nature of the gap and its anisotropy. Fig. shows the 
DFT±DMFT spectrum of the two most relevant orbitals, 
\mj\ = 1/2 and \mj\ = 3/2 of the Ce-J = 5/2 manifold. 
There is a substantial broadening due to the imaginary 
part of the self energy which is still huge at 290K, in¬ 
dicating that we are still considerably above the Kondo 
temperature. 

To get a better impression of the correlated bandstruc¬ 
ture, we hence set the imaginary part of the self energy to 
(essentially zero) in Fig.[^ This way one can also mimick 
the low temperature behavior since the imaginary part 
of the self energy vanishes like ^ around the Fermi 
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FIG. 5: Orbitally resolved DMFT self energies for real fre¬ 
quencies (jj at inverse temperatures /3 = 10 eV“^ and 40 eV“^. 
Top: real part; bottom: imaginary part. 



FIG. 6: DMFT spectra obtained by the maximum entropy 
method for T = 1160 K (middle) and 290 K (bottom) com¬ 
pared to DFT (top). 


level. One should keep in mind however, that the further 
build-up of the Kondo resoance might still reshuffle the 
bandstructure at lower temperatures. 



FIG. 7: DMFT spectra at T = 290 K resolved for the J = 5/2 
orbitals (left) and corresponding contributions from the other 
manifolds (right). 

Given the small crystal field splitting between |mj| = 
1/2 and |mj| =3/2 and its reversal by correlations as 
well as the much higher Kondo temperature for a SU(4) 
than a SU(2) Kondo effect, we can expect that both 
states will still contribute to the low temperature Kondo 
resonance. 

In principle, Fig.j^shows a Kondo insulator with a gap 
through most of the Brillouin zone. However, (i) along 
the 2 : direction from G to X there a crossing of bands and 
hence no gap. Also in the other directions, there is (ii) a 
negative indirect band gap. These are two independent 
meachisms of why CeRu 4 Sn 6 remains metallic with quite 
some anisotropy. 

The origin of the gap is the usual Kondo insulator sce¬ 
nario, i.e., the hybridization between Ce-/ and Ru-4d(Sn- 
5p) conduction bands. This can be seen by looking at the 
N and point in Fig. where there is a minimum of 
a (spin-split) conduction band around -0.2 eV of mainly 
Ru-4d character. Around the Fermi level these band go 
over smoothly into a flat predominantly Ce |mj| = 1/2 
band. In principle the same picture can also be found 
above the Fermi level. However, here it is somewhat less 
clear, since the conduction band maxima e.g. between 
iFg and G or iFs and G are not well separated from the 
other bands. 

Around X and iFs, the whole low-energy bandstruc¬ 
ture is shifted upwards. Here the gap is well above the 
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FIG. 8: DMFT /c-resolved spectra of the two most important 
orbitals, the J = bl2,mj = 3/2 (top) and J = 5/2, mj = 1/2 
(bottom) orbital, at T=290 K. White lines: DFT bandstruc- 
ture for comparison. 


Fermi level. Hence, we only have a direct band gap at 
most /c-points whereas the indirect gap is zero, we have 
a Kondo semimetal. Note that the term semimetal has 
been used before, see e.g. Refs.^^^K^, to classify CeNiSn 
but to the best of our knowledge in another connotation, 
i.e., to denote the weak metallic behavior - even with 
the Ikeda-Miyake scenaricPIin mind, i.e., a hybridization 
gap with nodes and hence a crossing of conduction and 
/ band. 

Moreover, along the 2 : direction, we even have a gap¬ 
less crossing of bands, see Fig. From Fig. |^we infer 
that at G the lowest (spin-split) band is \mj\ = 1/2, 
the second one \mj\ = 3/2 and the topmost the pre¬ 
dominately Ru conduction band. Hence, this crossing is 
between \mj\ = 1/2 and \mj\ = 3/2. 

Altogether this scenario is much more complicated 
than a simple model with flat /-band. In Ref. |4] it 
was concluded that since \mj\ = 1/2 orbitals have no 
nodes, the \mj\ = 1/2 Kondo insulator has to be gapped 


FIG. 9: DMFT /c-resolved spectrum (top) and that for an ar- 
tihcial S(ct;) = 5R(ST=290KDMFT(ct;)) + z0.005[eV] (bottom). 
White lines: DFT bandstructure for comparison. 


throughout the Brillouin zone. This scenario is however 
too simple for CeRu 4 Sn 6 . First of all, there are both, 
\mj\ = 1/2 and \mj\ =3/2, bands in the vicinity of the 
Fermi level; and we (i) hence observe a crossing of bands. 
Even if we had no such crossing, we would (ii) nontheless 
still have a Kondo semimetal with remnant metallicity 
instead of a Kondo insulator since the indirect band gap 
is negative. Because of these two reasons, the idealized 
modelling of Ref. |4] falls short of describing CeRu 4 Sn 6 . 


III. CONCLUSION 

In perspective of new x-ray absorption experiments^ 
we have presented further DFT-bDMFT results for 
CeRu 4 Sn 6 . In agreement with experiment, the dominat 
Ce orbital is J = 5/2, \mj\ = 1/2 which becomes the low¬ 
est crystal field level (and gets more occupied) when the 
Kondo effect sets in. However, there is also a secondary 
J = 5/2, \mj\ = 3/2 contribution to the Kondo reso- 

















6 


nance at the Fermi level. Because we have these two lev¬ 
els there is a band crossing along the z direction. Besides, 
the indirect Kondo gap is negative. That is even without 
band crossing we would still have a Kondo semimetal. 
Both mechanism on their own yield a remnant metallic- 
ity and anisotropy. 
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